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Abstract 

A full Navier-Stokes solver has been used to 
model transonic flow over three airfoil sections. 
The method uses a two-dimensional, implicit, con- 
servative finite difference scheme for solving the 
compressible Navier-Stokes equations. Results are 
presented as prescribed for the Viscous Transonic 
Airfoil Workshop to held at the AIAA 25th Aerospace 
Sciences Meeting. The NACA 0012, RAE 2822 and 
"Jones" airfoils have been investigated for both 
attached and separated transonic flows. Predic- 
tions for pressure distributions, loads, skin 
friction coefficients, boundary layer displacement 
thickness and velocity profiles are included and 
compared with experimental data when possible. 
Overall, the results are in good agreement with 
experimental data. 

Nomenclature 

a angle of attack, a (deg) 

c chordlength 

CD drag coefficient 

CF skin friction coefficient base on 


R average distance to outer boundary from 

leading edge 

RE Reynolds number based on chord 

s arclength of a grid line in the 

n-di recti on 

s e arclength of a grid line from the airfoil 

surface to the node with a velocity of Ue 

u,v cartesian velocities normalized by the 
speed of sound 

Ue maximum total velocity from the airfoil 

surface along a grid line 

Uinf free-stream velocity 

x,y cartesian coordinates normalized by 

chordlength 

XS shock location on the upper surface 

(chords) 

C chordwise direction of transformed 

coordinate system 



free-stream conditions 

n 

normal direction of transformed 
coordinate system 

:l 

lift coefficient 

p 

density 

CLmax 

maximum lift coefficient 

T 

time variable 

CP 

DCL/DA 

pressure coefficient 

slope of linear portion of CL versus a 

An 

distance of the first grid line off the 
airfoil surface (chords) 


DEL inverse of the number of upper surface 
grid points 

DSTAR boundary layer displacement thickness 
(chords ) 



e total energy of the fluid per unit volume 

J Jacobian of transformation 

Minf free-stream Mach number 

♦Graduate Research Assistant; student member AIAA. 
♦♦Associate Professor; member AIAA. 


Introduction 

In recent years, many computer programs have 
been written that predict the transonic flow prop- 
erties around an airfoil. The purpose of these 
codes is to help in the design of airfoils, such 
as those found on helicopters, propellers, wings, 
and turbomachinery. Accurate prediction of the 
flow field is necessary for performance, acoustic 
and aeroelastic analyses. In order to determine 
the quantitative ability of these codes, a workshop 
has been organized to compare a variety of flow 
solvers. These are all capable of producing air- 
foil solutions over a wide range of flow condi- 
tions, including both attached and separated 
transonic flow. The predictions from individual 
flow solvers are to be presented at the Viscous 
Transonic Airfoil Workshop in coordination with the 
AIAA 25th Aerospace Sciences Meeting. The results 
will be combined for the Fluid Dynamics Meeting six 
months later. The present study shows results from 
one of the flow solvers to be presented at the 
workshop. An unsteady, two-dimensional, full 
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Navier-Stokes solver is used to compute the flow- 
field around the three airfoils shown in Fig. 1. 
The required plots for the workshop are summarized 
in Table 1. Results for these test cases are 
presented to help determine the quantitative abil- 
ity of the present Navier-Stokes solver by com- 
paring the results with experimental data. 

Numerical Procedure 


E y F + C U G + c t q 


G = 


n v F + n G + n q 


R = 


K R + K S 
x y 


The program addresses the two-dimensional flow 
problem by solving the unsteady, two-dimensional, 
Reynolds-averaged, compressible Navier-Stokes 
equations on a body-fitted coordinate system in 
strong conservation form using an ADI procedure. 

The formulation has been described in Ref. 1, and 
only a brief outline is given here. All of the 
calculations were performed in a body-fitted coor- 
dinate system (s,o,t) which is linked to the 
cartesian coordinates according to the following 
one-to-one relationships: 

E = c(x,y,t) 
n = n(x,y,t) 

x = t (1) 

The Jacobian of the transformation is given by 

J = E x n y - n x S y = _ x ^) (2) 

and the metrics of the transformation are given by 
the relationship: 


Ey = -Jx„ 


Standard central differences were used to compute 
the quantities such as x^, y 5 etc., and 
these quantities in turn were used in Eqs. (2) and 
(3) to compute c x , e y , etc. At the 
boundaries, three-point one-sided differences were 
used to compute the metrics. The far field bound- 
aries are assumed to be undisturbed; except at the 
downstream boundary, where the velocities (u,v) 
and the entropy are extrapolated from the interior. 


In this coordinate system, the two-dimensional 
unsteady Navier-Stokes equations may be written as 
follows: 


where 

. A A A A 

dt + F 5 + G n = R c + S n 

(4) 


q = J~l(p, pu, pv, el 

(5) 

and p 

is the fluid density; u and 

v are the 


cartesian components of fluid velocity; e is the 
total energy^ of^thg fluid per unit volume. The 
quantities F, G, R, and S are given by: 


S = 


E R + E„S 


( 6 ) 


The terms F, R, etc. are the standard flux 
and viscous stress terms along the x-direction, in 
the original cartesian form. The terms G and 
S are likewise the flux and viscous terms along 
the y-direction. 

An eddy viscosity concept is used to model the 
turbulent momentum and energy transfer. The entire 
flow field is assumed to be turbulent. 

Since Eqs. (4) and (5) are coupled to each 
other, and are highly nonlinear, a stable, effi- 
cient solution procedure is required. In the 
present work the Beam-Warming algorithm, 5 as 
implemented by Steger,° was followed with some 
changes made to the artifical dissipation terms to 
avoid shock overshoots, and to avoid excessive 
smearing of embedded vortices. The viscous terms 
were lagged by one time step, while other terms in 
the governing equations were treated implicity, at 
the current time level. The nonlinear terms are 
split, with parts treated either implicitly or 
lagged in time. The details of the discretization 
procedure, the artifical dissipation terms used and 
several benchmark calculations are presented in 
Ref. 1. 


The original algebraic grid generator used by 
Sankar and Tang* has been replaced with the GRAPE 
code (GRids about Airfoils using Poisson's Equa- 
tions) developed by Sorenson. 2 This qives the 
Navier-Stokes solver greater flexibility to inves- 
tigate the effects of grid changes, as required by 
the workshop. For most cases, a 157 by 58 C-grid 
is used with 60 nodes defining the wake. This 
leaves 97 nodes wrapped around the upper and lower 
surfaces. The distance off the airfoil surface to 
the first grid line is 0.00005 chordlengths. The 
GRAPE code allows the user to vary the clustering 
of the grid lines near the airfoil surface. A 
typical base grid used in this study is shown in 
Fig. 2. 

All runs were done on the CRAY-XMP at the NASA 
Lewis Research Center. The number of iterations 
required for convergence range between 1500 and 
5000, depending on the flow conditions. The res- 
pective vectorized CPU times were 519 and 1730 sec. 
The calculations continue until the spatial- 
maximum density residual drops by four orders of 
magnitude from the initial condition. The solver 
has been vectorized so that the ratio of scalar 
CPU time to vectorized CPU time is about 2.68. The 
total calculation time is approximately 3.8x10“^ 
sec per point per iteration. The maximum memory 
required is 960 kwords using a 157 by 58 grid. 
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Results 

Case ft: NACA 0012 Airfoil 

A study of the grid is necessary to determine 
the sensitivity of the predictions to variations 
in the grid. The present results include both grid 
refinement and boundary position comparisons. 

Define a grid refinement parameter, DEL, to be the 
inverse of the number of upper surface grid points. 
For the 157 by 58 base grid, the value of DEL is 
0.021. A variation in DEL represents a propor- 
tional change in the grid size from the base grid. 
The outer boundary position remains fixed at six 
chordlengths. (This ratio is equally enforced in 
the wake region.) Figures 3 and 4 show lift and 
drag predictions as a function of DEL for the 
test case: MINF = 0.70, a = 1.49° and RE = 9.0 
million. The asymptotic behavior for loads as DEL 
becomes small suggests that the 157 by 58 base grid 
gives reasonable predictions. 

Next, the outer boundary position (R) was 
varied from 2 to 20 chordlengths away from the 
leading edge using a 157 by 58 grid. Figures 5 and 
6 justify the base grid results for lift and drag 
predictions. When R equals 15 and 20, the grid 
spacing near the airfoil surface increases signi- 
ficantly. This explains the nonasymptotic behavior 
of the drag predictions in Fig. 6 as R increases. 

The distance of the first grid line off the 
airfoil surface (an) must be carefully chosen to 
properly model the boundary layer. Figure 7 shows 
how this parameter can effect the pressure distri- 
butions for three values of an: 0.00002, 

0.00005, and 0.0005. When an = 0.0005, the 
shock is smeared and the magnitude of pressure is 
higher aft of the shock. A value of an = 0.00005 
appears to give satisfactory predictions and is 
used in the base grid. The base grid (157 by 58, 

R = 6 chords, an = 0.00005 chords) is now defined 
for the NACA 0012 airfoil and will be used for the 
remainder of the test cases. 

Figures 8 through 10 show pressure distribu- 
tions over the NACA 0012 airfoil for several ranges 
of velocity and angle of attack. Results for two 
values of angle of attack are plotted along with 
experimental data from Harris. ^ The higher angle 
corresponds to the value used in the experiment, 
where the lower angle is corrected for wind tunnel 
wall interference (as suggested in Ref. 3). 

Results show that wind tunnel wall interference can 
significantly change the behavior of the pressure 
distributions. The corrected angle of attack tends 
to give better agreement with experiment. The 
overall agreement is very good for test cases in 
Figs. 8 and 9. In the case of a strong shock, the 
predicted shock location is about 0.10 chordlengths 
aft of the experimentally determined shock location 
(see Fig. 10). The flow solver requires as input 
a value for the artificial (numerical) damping. 

It is important to minimize this value in order to 
properly model the Navier-Stokes equations. For 
the test case in Fig. 10, higher numerical damping 
tends to make the predicted shock location appear 
to be closer to the experimental shock location. 
Higher numerical damping also results in more shock 
smearing. 

A plot of CL versus o (a) is presented in 
Fig. 11 for Minf = 0.70 and RE = 9.0 million. 


The predictions are good for flow conditions where 
experimental data exists. A variable time step is 
used in steady-state solutions to decrease the 
convergence time. However, when unsteady flow 
exists, a time accurate solution is required. Two 
time accurate predictions are included in Fig. 11 
for a equal to 5.86 and 6.5°. Unsteady flow due 
to separation is expected at these angles, but 
predictions between the variable time step and time 
accurate cases appear to be nearly identical. 

Plots of CL versus CD and CD versus Minf are 
shown in Figs. 12 and 13. The overall agreement 
with experiment is good. 

Figure 14 shows how the shock location moves 
aft as the free-stream Mach number increases 
(a = 0.0°, RE = 9.0 million). Predictions for 
CLmax versus Minf are presented in Fig. 15. All 
runs used to determine CLmax were time accurate 
with increments for a of 0.5°. The magnitude of 
the prediction for CLmax is expected to be some- 
what inaccurate due to this discretization of a. 
Figure 16 shows how the slope of the linear portion 
of CL versus a is predicted to change with 
Minf. 

Case B: RAE 2822 Airfoil 

Results from two test cases are presented for 
the RAE 2822 airfoil and compared with experimental 
data from Cook, McDonald and Firman. 4 The effect 
of the wind tunnel interference is not easily 
determined. Hence, a range of flow conditions were 
tried to match the pressure distributions. A 
"matched" pressure distribution is the best agree- 
ment of the predictions with the experimental data 
that can be found by adjusting the Mach number and 
the angle of attack. In cases where a strong shock 
exists, a good "matched" pressure distribution is 
difficult and only an overall agreement can be 
estimated. The experimental pressure distributions 
are shown in Fig. 17 for Minf = 0.725, a = 2.92° 
and RE = 6.5 million. Predictions for 
Minf = 0.73 and a = 2.79° are plotted from the 
flow solver using a 157 by 58 base grid similiar 
to the base grid used in the NACA 0012 airfoil 
analysis. (The conclusions from the grid study 
done for the NACA 0012 airfoil were assumed to be 
the same for both the RAE 2822 and Jones airfoils). 
The agreement is very good for most regions on the 
airfoil. The expansion near the leading edge on 
the upper surface shown in the experimental data 
is probably due to the transition location being 
fixed at x/c = 0.03. This is not modeled in the 
flow solver since the entire flow field is assumed 
to be turbulent. The integrated lift and drag 
coefficients from the experiment are respectively 
0.7433 and 0.0127, which correspond to 0.7423 and 
0.0134 from the Navier-Stokes code. 

The calculation of displacement thickness, 
DSTAR, is not straight forward. The definition of 
edge velocity, Ue, is not as well defined for an 
airfoil as it is for a flat plate. Separated flow 
and curvature from the inviscid velocity variation 
further complicates the computation. For the 
present analysis, the edge velocity used to cal- 
culate DSTAR is the first maximum total velocity 
away from the airfoil surface. The boundary layer 
is the arclength of the grid line measured from the 
airfoil surface to the node identified as having 
the velocity, Ue. Figure 18 shows the predicted 
DSTAR along the upper surface of the current test 
case. The boundary layer becomes thicker aft of 
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the shock, as expected. The skin friction coeffi- 
cient decreases along the upper surfaces, as shown 
in Fig. 19. Velocity profiles are shown for two 
locations on the upper surface: x/c = 0.319 and 

x/c = 0.956. The magnitudes of the velocities used 
to calculate U/Uinf are taken from the component 
of velocity in the free-stream direction. The 
results are presented in Figs. 20 and 21. 


The second test case for the RAE 2822 airfoil 
corresponds to the following experimental flow 
conditions: MINF = 0.75, a = 3.19° and 

RE = 6.2 million. The "matched" conditions 
selected^for the flow solver are: Minf = 0.755 and 

a = 3.03°. A comparison of the pressure distribu- 
tions is shown in Fig. 22. Again, minor discrep- 
ancies exist for the upper surface pressure near 
the leading edge. The shock location is predicted 
slightly aft of the experiment. Recall that this 
trend was also observed for NACA 0012 airfoil for 
the test case with a strong shock (see Fig. 10). 

The flow behind a strong shock is unsteady due to 
separation, which challenges the turbulence model 
to properly estimate the eddy viscosity. Coakley^ 
has shown how important the turbulence model can 
be for this test case. He found that predictions 
using the Baldwin-Lomax model (used in the present 
solver) causes the shock to be located aft of 
experiment, while predictions using other turbul- 
ence models improve the agreement. 


The displacement thickness of the boundary 
layer on the upper surface is shown in Fig. 23. 
Agreement between analysis and experiment is good 
in the region forward of the shock. Figure 24 
shows the corresponding skin friction coefficient 
distribution. Two separated regions are predicted 
near the shock and near the trailing edge. 

Velocity profiles at x/c = 0.75 and x/c = 0.897 
are presented in Figs. 25 and 26. The predictions 
appears to show attached flow, but the scale 
chosen for the plots does not show that the flow 
is predicted to be separated very close to the 
airfoil surface. The separation regions near the 
airfoil surface are shown best in the skin fric- 
tion coefficient distribution. The sonic line 
prediction is shown if Fig. 27 as a dotted line. 

Case C: Jones Airfoil 


The last airfoil investigated in this study 
is the Jones "supercritical-type" airfoil, devel- 
oped by R.T. Jones of the NASA Ames Research 
Center. This airfoil was selected for the Viscous 
Transonic Airfoil Workshop because it has not been 
widely tested either numerically or experimentally 
and serves as a good "blind" calculation for flow 
solvers. 

Predictions for CL versus a and CL 
versus CD are shown in Figs. 28 and 29. The 
lift curve slope begins to drop, but never drops 
enough to give a good CLmax prediction. The 
solution residuals become oscillatory at higher 
angles of attack, which usually indicates unsteady 
flow. Time accurate solutions were tried for 
a = 6.5 and a = 7.0°. A slight drop in the lift 
coefficients is observed, as shown in Fig. 28. A 
plot of predicted CD versus Minf is shown in 
Fig. 30 for a = 0.0°. 

A supercritical test case was run for 
Minf = 0.75, a = 2.0° and RE = 9.0 million. The 
pressure distributions in Fig. 31 shows a strong 


shock predicted on the upper surface near the 
forty-percent chord location. The boundary layer 
thickens near the shock, as shown in Fig. 32. The 
sudden rise in the displacement thickness near the 
trailing edge is due to the rapid change in the 
edge velocity location in the chordwise direction. 
The method for calculating displacement thickness 
uses the first maximum total velocity off the air- 
foil surface. If the velocity profile is very 
complicated, as it is near the trailing edge, the 
velocity is not easily determined. The plot is 
dotted in this region to represent the uncertainty 
of the DSTAR calculation. Skin friction coeffi- 
cients and velocity profiles at x/c = 0.613 and 
x/c = 0.897 are presented in Figs. 33 through 35. 

Results are also presented for a separated, 
supercritical case (Minf = 0.85, a = 2.0°) and 
summarized in Figs. 36 through 40. In this case, 
shocks are predicted on both the upper and lower 
surfaces. Separation occurs from about x/c = 0.80 
to the trailing edge on the upper surface. Again, 
the displacement thickness calculations in this 
region are not well defined. 

Conclusion 


A full Navier-Stokes flow solver has been 
evaluated for modeling transonic flow over the 
NACA 0012, RAE 2822 and Jones airfoil sections. A 
157 by 58 C-grid, with an average outer boundary 
distance of six chords appears to be adequate for 
most flow conditions. Comparisons with experi- 
mental data have been shown for both the NACA 0012 
and RAE 2822 airfoils. In cases with strong shocks 
and separated flow, the predicted shock locations 
are found to be aft of the experimentally deter- 
mined locations. Separation near the shock and the 
trailing edge is also predicted for these cases. 
Results for lift coefficient versus a predict a 
loss in lift at higher angles of attack and require 
a time-accurate solution when the flow is unsteady. 
Wind tunnel interference corrections show improved 
agreement between the predicted and experimental 
pressure distributions. Overall, the results are 
in good agreement with experimental data. 
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TABLE 1. - TEST CASES SPECIFIED BY THE VISCOUS 
TRANSONIC AIRFOIL WORKSHOP 


NACA 0012 airfoil (RE = 9.0 million) 

CP versus 

x/c 

Minf = 0.70, a = 1.86° 

CP versus 

x/c 

Minf = 0.55, a = 9.86° 

CP versus 

x/c 

Minf = 0.799, a = 2.86° 

CL versus 

a 

Minf = 0.70, a = 1,3, 5, 6, 7° 

CL versus 

CD 

Minf = 0.70, a =0,1, 3, 5, 6, 7° 

CD versus Minf 

a = 0, Minf = 0.5 to 1.1 

CL versus 

DEL 

Minf = 0.70, a = 1.86° 

CD versus 

DEL 

Minf = 0.70, a = 1.86° 

CL versus 

R 

Minf = 0.70, a = 1.86° 

CD versus 

R 

Minf = 0.70, a = 1.86° 

XS versus 

Minf 

a = 0, Minf = 0.8 to 1.1 

CLmax versus Minf 

Minf = 0.3, 0.5, 0.6, 0.7, 0.8 

DCL/DA versus Minf 

a = 0, Minf = 0.5 to 1.1 


RAE 

2822 Airfoil 

CP versus 

x/c 

Minf = 0.725, a = 2.92° 



RE = 6.5 million. 

DSTAR/c versus x/c 

upper surface 

CF versus 

x/c 

upper surface 

y/c versus 

U/Uinf 

x/c = 0.319, upper surface 

y/c versus 

U/Uinf 

x/c = 0.95, upper surface 

CP versus 

x/c 

Minf = 0.75, a = 3.19 , 



RE = 6.2 million 

DSTAR/c versus x/c 

upper surface 

CF versus 

x/c 

upper surface 

y/c versus 

U/Uinf 

x/c = 0.75, upper surface 

y/c versus 

U/Uinf 

x/c = 0.90, upper surface 

Sonic line 




Jones airfoil (RE = 9.0 million) 


CL 

versus 

a 

Minf = 0.7, a = 0,1, 3, 5, 6, 7° 

CL 

versus 

CD 

Minf = 0.7, a = 0,1, 3, 5, 6, 7 

CD 

versus 

Minf 

a = 0, Minf = 0.5 to 1.1 

CP 

versus 

x/c 

Minf = 0.75, a = 2.0° 

DSTAR/c versus x/c 

upper surface 

CF 

versus 

x/c 

upper surface 

y/c 

versus 

U/Uinf 

x/c = 0.6, upper surface 

y/c 

versus 

U/Uinf 

x/c = 0.9, upper surface 

CP 

versus 

x/c 

Minf = 0.85, a = 2.0° 

DSTAR/c versus x/c 

upper surface 

CF 

versus 

x/c 

upper surface 

y/c 

versus 

U/Uinf 

x/c = 0.6, upper surface 

y/c 

versus 

U/Uinf 

x/c = 0.9, upper surface 
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FIGURE 1.- AIRFOIL SECTIONS. 



FIGURE 2. - 157X58 BASE GRID USED IN COMPUTATIONS, AVERAGE OUTER BOUNDARY DISTANCE 
EQUALS SIX CHORDS. 
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FIGURE 7.- EFFECT OF GRID SPACING NORMAL TO THE 
AIRFOIL SURFACE ON PRESSURE DISTRIBUTIONS. NACA 
0012 AIRFOIL, MINF = 8.80, Re = 9.00X10 6 . 



FIGURE 8. - PRESSURE DISTRIBUTIONS ON 
NACA 0012 AIRFOIL. MINF = 0.70; 

A = 1.86 (1.A9); Re = 9.00x10. 
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FIGURE 14. - PREDICTED SHOCK WAVE LOCATION 
VERSUS MINF. NACA 0012 AIRFOIL; A = 0; 
Re = 9.00X10 6 . 
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FIGURE 19. - SKIN FRICTION COEFFICIENT ON 
UPPER SURFACE OF RAE 2822 AIRFOIL. 

A = 2.92; MINF = 0.725; Re = 6.50x10 6 . 



U/UINF 

FIGURE 20. - VELOCITY PROFILE AT 
X/C = 0.319. UPPER SURFACE OF 
RAE 2822 AIRFOIL; A = 2.92; 
MINF = 0.725; Re = 6.50X10 6 . 
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FIGURE 21. - VELOCITY PROFILE AT 
X/C =0.95. UPPER SURFACE OF 
RAE 2822 AIRFOIL; A = 2.92; 
MINF = 0.725. Re = 6.50X10 6 . 
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FIGURE 22. - PRESSURE DISTRIBUTIONS ON RAE 

2822 AIRFOIL. MINF = 0.750; A = 3.19; 

Re = 6.20X10 1 * * * * 6 . 
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FIGURE 23.- DISPLACEMENT THICKNESS ON UPPER 
SURFACE OF RAE 2822 AIRFOIL. A = 3.19, 
MINF = 0.750, Re = 6.20X10 6 . 



FIGURE 24. - SKIN FRICTION COEFFICIENT ON 
UPPER SURFACE OF RAE 2822 AIRFOIL. 

A = 3.19; MINF = 0.750; Re = 6.20x10 6 . 




FIGURE 25. - VELOCITY PROFILE AT 
X/C = 0.750. UPPER SURFACE OF 
RAE 2822 AIRFOIL; A * 3.19; 
HINF = 0.750; Re = S.20X10 6 . 



FIGURE 26. - VELOCITY PROFILE AT 
X/C = 0.90. UPPER SURFACE OF 
RAE 2822 AIRFOIL; A = 3.19; 
MINF = 0.750; Re = 6.20X10 6 . 




FIGURE 27. - SONIC LINE PREDICTION FOR RAE 2822 AIRFOIL. 
A = 3.19; MINF = 0.75; Re = 6.20X10 6 . 
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FIGURE 28. - PREDICTED CL VERSUS ALPHA. 
JONES AIRFOIL; MINF = 0.70; 

Re = 9.00X10 6 . 
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FIGURE 32. - PREDICTED DISPLACEMENT THICK- 
NESS ON UPPER SURFACE OF JONES AIRFOIL. 

A = 2.00; MINF = 0.75; Re = 9.00X10 6 . 
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FIGURE 33. - PREDICTED SKIN FRICTION 
COEFFICIENT ON UPPER SURFACE OF JONES 
AIRFOIL. A = 2.00; MINF = 0.75; 

Re = 9.00x10**. 
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FIGURE 36. - PREDICTED PRESSURE DISTRIBU 
TIONS ON JONES AIRFOIL. MINF = 0.85; 

A = 2.00; Re = 9.00x10 . 

















FIGURE HO. - PREDICTED VELOCITY PROFILE 
AT X/C = 0.9. UPPER SURFACE OF JONES 
AIRFOIL; A = 2.00; MINF = 0.85; 

Re = 9.00x10®. 


■ 


1. Report No. j(y| 88912 2. Government Accession No. 

AIAA-87-0420 

3. Recipient’s Catalog No. 

4. Title and Subtitle 

Analysis of Viscous Transonic Flow Over Airfoil 
Sections 

5. Report Date 

6. Performing Organization Code 

535-03-01 

7. Authors) 

Dennis L. Huff, Jlunn-Chl Wu, and L.N. Sankar 

8. Performing Organization Report No. 

E-3340 


10. Work Unit No. 

9. Performing Organization Name and Address 


National Aeronautics and Space Administration 
Lewis Research Center 

11. Contract or Grant No. 

Cleveland, Ohio 44135 

13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

— 

Technical Memorandum 

14. Sponsoring Agency Code 


15. Supplementary Notes 


Prepared for the 25th Aerospace Sciences Meeting, sponsored by the American 
Institute of Aeronautics and Astronautics, Reno, Nevada, January 12-15, 1987. 
Dennis L. Huff, NASA Lewis Research Center; Jlunn-Chl Wu and L.N. Sankar, 

Georgia Institute of Technology, School of Aerospace Engineering, Atlanta, 

Georgia 30332. 

16. Abstract 

A full Navler-Stokes solver has been used to model transonic flow over three 
airfoil sections. The method uses a two-dimensional. Implicit, conservative 
finite difference scheme for solving the compressible Navler-Stokes equations. 
Results are presented as prescribed for the Viscous Transonic Airfoil Workshop 
to be held at the AIAA 25th Aerospace Sciences Meeting. The NACA 0012, RAE 2822 
and "Jones" airfoils have been Investigated for both attached and separated 
transonic flows. Predictions for pressure distributions, loads, skin friction 
coefficients, boundary layer displacement thickness and velocity profiles are 
Included and compared with experimental data when possible. Overall, the results 
are In good agreement with experimental data. 


17. Key Words (Suggested by Authors)) 18. Distribution Statement 

Navler-Stokes equations; Transonic flow; Unclassified - unlimited 
Airfoils; Computational fluid dynamics STAR Category 34 


19. Security Classlf. (of this report) 20. Security Classlf. (of this page) 21. No. of pages 22. Price* 

Unclassified Unclassified 


For sale by the National Technical Information Service. Springfield, Virginia 22161 


















